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Abstract. We present a novel quantum algorithm for classification of 
images. The algorithm is constructed using principal component anal¬ 
ysis and von Neuman quantum measurements. In order to apply the 
algorithm we present a new quantum representation of grayscale images. 


1 Introduction 


At the end of the last century a new paradigm of computation was proposed i.e. 
quantum computation. Although it is not yet obvious whether useful quantum 
computers can be constructed, the field of quantum algorithms development in 
recent years progresses very rapidly [2], [1]. For example many new algorithms for 
quantum machine learning and quantum image processing were recently created 
[9], [11]. 

In this work we introduce an algorithm for image classification of grayscale 
images based on classical principal component analysis (PCA) and quantum 
measurement. The general idea behind the algorithm is following. Given a set 
of training images, using PCA we train a classifier to detect images similar 
to those in the training set. Effectively we divide the image signal space into 
two orthogonal subspaces. The first one — spanned by the leading principal 
components — catches the most of the variability of the signal in the training 
set, the second one consists mostly of noise. 

After the classifier is constructed the leading principal components are used 
to create a projector onto a subspace of quantum states. The image which is 
being classified is also encoded on a quantum state, and then measured using 
the projector defined above. 

The paper is organised as follows in Section 2 we recall basic notions of 
quantum computation, in Section 3 we shortly discuss state of the art in quantum 
image processing, in Section 4 we introduce image classification algorithm. And 
finally in Section 5 we draw conclusions. 


2 Essentials of quantum computation 


Lets consider the most basic model of a quantum system - a qubit - elementary 
quantum system with two basic physical states. In order to provide mathemat¬ 
ical description of a state of a qubit we choose an orthonormal basis in the 
corresponding Hilbert Space. In this case we consider two dimensional Hilbert 
Space. Our basis will consist two vectors that in the braket notation take the 
form 
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The \x) vector is called ‘kef and its Hermitian conjugation (|a;))'^ = {x\ is called 
‘bra’. We can represent any valid state of a qubit \'>p) as normalized linear com¬ 
bination of the basis vectors: 

|V’) = ai|0)-|- \-ak\k), (2) 
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where ai, ..., G C and ^ = 1. 
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The operation which allows us to join k independent qunit systems is the 
tensor product. Lets take k qunit states 
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We can write their joint state in as 





We can also consider more general quantum systems called qudits. Let ji/') 
(ket) be a normed column vector from Hilbert space C" with orthonormal basis 
{|*)}?=i- Dual vector to ket is (f/'l (bra). In such case the state of the system is 
represented as If/') = V'il*)- 

n 

We denote inner product of jf/') and \(j)) by (^|f/') = X) It has three 

properties: 

1 . (f/’|f/’) > 0 where equality holds iff jf/') = 0 , 

2. ((/)|f/^) = (f/'|</))*, 

3. (f/>|(ai|(/)i) + a2|(/>2)) = ai(f/'|</'i) + a2{ip\(j)2)- 

n n 

Furthermore |f/’)(^f'| = X X l*)(j| ’"'ih be their outer product. 

i=l j=l 

One of the most important concepts in quantum information is the mea¬ 
surement. The mathematical model of a measurement is as follows. At first we 
define a set of outcomes F. Then we assign corresponding measurement opera¬ 
tors We request that the measurement operators satisfy the condition 

= Pj and X 7 P '1 = I- 

The probability that we obtain outcome 7 when measuring a state \(j)) is 
equal to 


Prh,\(l))) = {(l)\P'i\4')- ( 6 ) 

If we instantly measure the system for the second time, the outcome willstill be 
equal to 7 with certainty because after the first measurement the state of the 
system changes into a state 

{<p\p,W/^- 

3 Quantum image processing — state of the art 

There are various ways in which classical data can be encoded on quantum states. 
The specific encoding depends on the type of the data and quantum algorithms 
that one wishes to execute. 

3.1 Quantum representations of digital images 

Below we recall various representations of quantum images proposed in recent 
years. 

In the Qubit Lattice representation of grayscale images proposed in [14] the 
intensity of pixel at position y,x is encoded on qubit \q)y,x- 

The Real Ket representation introduced in [ 6 ] stores 2" x 2" grayscale images 
in unnormalised quantum states of the form 
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where G R and subsequent qu4its serve as the position of pixel encoded 

in a quad-tree. 

Flexible representation of quantum images (FRQI) captures information about 
pixel colours and their corresponding position. It is inspired by the pixel repre¬ 
sentation for images in the classical computers. The information is gather into 
a quantum state dehned as follows 
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where 9i G [O, and constitutes the vector encoding colours, |0), |1) is a fixed 
basis of a two dimensional complex Hilbert space, and |z) is a basis of 2^” di¬ 
mensional space responsible for encoding position in the image. The colour is 
encoded in a 2D vector by cos di |0) -fsindi |1) which is connected by a Kronecker 
product with a vector |z) responsible for a position in the image. 

A novel enhanced quantum representation (NEQR) of digital images pro¬ 
posed in [18] encodes a grayscale 2" x 2” image in a quantum state of the form 

X—0 i—0 

where Cy^ is a discrete value of image intensity, quantised with q levels of quan¬ 
tisations of pixel at position {y,x). 

3.2 Quantum image processing algorithms 

The sub-field of quantum computation that deals with algorithms development 
for quantum image processing is developing very rapidly. At least a hundred 
papers discussing this subject were published in recent fifteen years. 

It should be noted that, some of classical image transformations already 
have their quantum analogues. For example we can list here quantum Fourier 
transform [10], quantum discrete cosine transform [5,12], and quantum wavelet 
transform [4]. 

There exists several clever techniques to process images encoded in quantum 
states for example in [3] authors propose a way to perform template matching 
algorithm using quantum Fourier transform and amplitude amplification. In pa¬ 
per [8] the authors extended the use of quantum circuit models for quantum 
image representation and processing. They developed three strategies to extend 
the number of geometric transformations [7] on quantum images using the FRQI 
representation of quantum images. In [16] authors propose quantum algorithms 
for edge detection and image filtering based on projective measurement. In [17] 
the authors propose a model for storing and operating on infra-red images. 

Complex quantum image processing requires a number of basic algorithmic 
primitives. In [15] authors develop quantum image translation, which maps the 
position of each picture element into a new position. In [19] an algorithm for 
comparing colour quantum images based on FRQI model is described. 



4 Algorithm 


The aim of our algorithm is classification of quantum images. The input of 
the algorithm is a quantum representation of an image which we want to test. 
Algorithm requires a set of principal components. The output is “yes” or “no” 
and answers the question whether the image exhibits features represented by 
principal components. 

4.1 Principal Component Analysis 

In order to create our quantum classifier we use Principal Component Analysis 
(PCA). This technique has been successfully applied in the domain of signal 
processing to various datasets. In celebrated classical paper [13] it was applied 
to classification of human faces. 

Suppose we have matrix of data A € Mm,n with rank k < m. The matrix is 
composed of vertically stacked horizontal sample vectors. We assume that our 
samples are normalised i.e. have I 2 norm equal to one. 

Then by SVD we have A = USV'^. where U G Mm and V € Mn are 
orthogonal matrices. The matrix S = diagjcri,... ,< 7 ^} is such that ui > a 2 > 

■ ■ ■> <Jk > cfk+i = ■ ■ ■ = CTg = 0, with q = min(m, n). 

The numbers ai are called singular values, i.e. non-negative square roots of 
the eigenvalues of AA"^. The columns of U are eigenvectors of AA^ and the 
columns of V are eigenvectors of A^A. The i-th column vector of the matrix U:^i 
is called the i-th principal component of the data. 

4.2 Quantum image representation 

Suppose we have a features vector of n values |A) = where Xi G [0,1]. 

The quantum system encoding the data from the feature space will be a direct 
sum H = (C*)®”. Quantum representation |^(A)) G "H of a picture |A) will be 
a mapping from [ 0 , Ij" to H defined by 
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|<i>(A)) = ^ 0 |</)(x.)), ( 8 ) 

where pixels are represented by 

\Hxi)) =x,\Q) + ^l-xl\l). (9) 

We will use quantum representation of vectors from PCA in the same way. 
Let {|Vi)}f=ibe a set of principal components with values vi^i G [—1,1] and 
{|^(Vi))}f^i be a set of quantum representations of them. Representation of 
each quantum vector chosen from PCA is encoded on different 2-dimensional 
subspace of C* such that each of the pixels is transformed into 

\<Pivj,i)) = w^,^|0) -f + 1), 


( 10 ) 



where j G {1,2,3,..., s} and i is a pixel index. The whole principal compo¬ 
nent representation is composed of the pixel representations the same way as 
in eq. (8). Lets take two vectors \Vj), \Vi) and their quantum representation 
l^(V^)), \^{Vi)). Inner product of \Vi) and \Vj) is 

n 

{Vi\V,) ( 11 ) 
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and for corresponding quantum representations one reads 
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where the first equality is from eq. (8) and the last one is implied by orthonor¬ 
mality of the basis vectors. From eq. (11) and eq. (12) we derive that inner 
product of two vectors is equal to inner product of quantum representations of 
these vectors with respect to a constant factor 1/n. This is an important feature 
of the introduced representation, which is significant for the algorithm. 


4.3 Construction of measurement 

The quantum algorithm for principal component analysis is based on classical 
methods for determining the characteristic subspace of the data set in the fea¬ 
tures space. Thus we take s principal components {|V;)}f^;^ that describe the 
data set crucial properties. The quantum algorithm for principal components 
analysis will utilise the system % dehned in the previous section. In order to 
use the classically computed components in the quantum algorithm we need to 
convert our principal components into quantum representation {|^(V/))})*^j^. 

The developed algorithm is based on the quantum measurement schema. We 
consider two elements output set F = {yes, no}. The first of the resulting labels 
will correspond to the principal components subspace and the other to the rest 
of the feature space. Thus we create two measurement operators 11 and 1 — 77. 
The principal components projection operator 77 is of the form 
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n = Y^\i>{Vi)){m)\- 
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4.4 Measurement probabilities 


Let \X) be a input feature vector and set of the principal components. 

In the classical model we measure the likelihood of the input data being in the 
control set in the following way 


M = Y^\{X\Vi)\\ (14) 

/=! 

Now let |^(-^)) be a quantum representation of the input and be 

quantum representation of principal components with projector 7T constructed 
as in eq. (13). Then the probability of the result of the measurement being Y 
for a given input is 


Pr(yes|X) = (<l>(X)|i7|<Z>(X)) = 1 ^ |(W|<I>(1^))P = (15) 

1=1 

where the last equation results from eq. (11) and eq. (12). Thus the probability 
Pr(yes|X) is linearly dependent on the classical likelihood measure M with 
respect to a factor 1/n^. where last equation is from eq. (11) and eq. (12). Thus 
the probability Pr(yes|X) linearly dependent on the classical likelihood measure 
M with factor 1/n^. 

Because of the factor l/n? we perform n? tests. We assume that we have 
copies of the quantum representation of the vector \X). We perform the 
measurement 17 on each of the copies. If any of the measurements returns “yes” 
then the algorithm returns positive answer. If not, the answer is negative. The 
probability that our algorithm will return the output “no” for a given input vector 
|X) is equal to 

Pr,n<iio\X) = (1 - Pr{jio\X)f. (16) 

Probability of positively classifying the input image in most of the cases is close 
to the classical likelihood measure. In general the probability is slightly lower. 
Thus the algorithm trifle favors the negative answer. 


5 Concluding remarks 

In this paper we provided a new quantum representation of digital images and 
an algorithm for classification of said images. The principal component analysis 
is used during the learning phase of the algorithm which is performed classi¬ 
cally and its goal is to construct a quantum measuring device. Classification is 
performed by applying the measurement apparatus on the quantum states that 
represent input images. The measurement is performed on multiple copies of the 
image. Therefore the paper provides complete system for classification of digital 
images. 
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